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ABSTRACT 

The  task  of  protecting  soldier  eyesight  and 
battlefield  optical  sensors  from  potentially  damaging 
effects  of  high-intensity  light  has  stimulated  interest  in 
electro-optic  (EO)  materials  such  as  EO  chromophores 
inserted  into  polymeric  material.  A  priori  design  of 
optimal  EO  chromophores  of  Army  relevance  requires 
the  capability  to  calculate  accurately  both  the  linear 
absorption  properties  (in  effect,  the  color)  and  non-linear 
optical  (NLO)  properties  of  EO  chromophores.  Here  we 
discuss  two  major  challenges  for  the  computational 
methods  used  in  EO  materials  research:  1)  efficiency  in 
exploring  an  extremely  large  space  of  possible  EO- 
molecular  architectures;  and  2)  accuracy  of  calculating 
ultraviolet-visible  (UV-Vis)  spectra.  We  show  that  the 
constrained  inverse  molecular  design  algorithm  (cIMD) 
and  the  long-range  corrected  density  functional  method 
(LC-DFT)  facilitate  the  discovery  of  new  chromophores 
for  Army  applications. 

1.  INTRODUCTION 

The  development  and  use  of  advanced  materials  in 
future  military  systems  are  critical  to  attaining 
performance  goals  of  the  Future  Force.  It  is  equally 
important  to  establish  safeguards  against  pervasive 
threats  by  developing  protective  EO  materials.  These 
materials  should  be  transparent  in  the  visible  and  possess 
NLO  properties  characterized  by  a  high  value  of 
molecular  hyperpolarizability  ( (3 )  (Pritchett  et  al.,2008). 


optimal  UV-Vis  absorption  spectra  of  chromophores 
with  large  (3  would  be  transparent  in  the  500  to  600  nm 
range  of  wavelengths  corresponding  to  a  yellow  spectral 
range.  This  requirement  can  be  achieved  if  the  main 
absorption  peak  Xmax»  700nm  and  the  potential  second, 
smaller  peak  is  at  Xmax-i  <  420nm. 

In  this  paper,  we  investigate  two  classes  of 
chromophores:  1)  tolane-  and  2)  FTC-based  (Fig.  1).  The 
tolane-based  chromophores  are  transparent  in  the 
desirable,  yellow  visible  range  but  have  only  mid-range 
p-values.  On  the  other  hand,  the  FTC-based  dyes  have 
very  large  P-values,  but  absorb  strongly  in  the  500-600 
nm  range.  Figure  1  shows  schematically  the  UV-Vis 
absorption  spectra  of  these  dyes  in  the  visible  range.  The 
challenge  is  to  modify  the  architecture  of  tolane  dyes  to 
increase  their  p-values  without  affecting  their  excellent 
transparency;  while  for  FTC  dyes  the  primary  goal  is  to 
enlarge  the  visibility  window  in  the  yellow  range. 


This  represents  a  considerable  development 
dilemma,  as  transparent  chromophores  are  typically 
small  molecules  with  low  EO  response.  Molecules  with 
the  desired  large  P-values  are  frequently  opaque  or,  at 
best,  have  small  windows  of  visibility,  and  are,  therefore, 
useless  for  Army  applications.  For  human  eyesight,  the 


Figure  1.  Structure  of  FTC  and  tolane  and  their 
absorption  spectra 

The  chromophore  with  phenyl- ethynyl  (tolane) 
motif  is  presented  in  Fig.  1.  A  systematic  study  of 
bis(phenylethynyl)  benzenes  with  various  donor  and 
acceptor  substitutions  has  been  reported  (Nguyen  et  al., 
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1997).  All  the  compounds  have  ^max  <  400  nm  and  the 
nitro-  aminomethyl  substitutions  possess  the  largest  (3  of 
57xl030  esu.  The  2-dicyanomethylen-3-cuano-4-{2-[E- 
(4-N,N-di(2-acetoxyethyl)-amino)-phenylene-(3 ,4- 
dibutyl)  thien-5]-E-vinyl}  -5, 5 -dimethyl-2, 5 -dihydrofuran 
(denoted  FTC)  chromophore  were  studied  extensively, 
both  experimentally  (Robinson  et  ah,  1999)  and 
theoretically  (Kinnibrugh  et  ah,  2006). 

Density  functional  theory  (DFT)  is  now  becoming 
the  most  popular  method  in  electronic-structure 
calculations  as  it  represents  a  good  compromise  between 
computational  performance  and  accuracy.  Recently,  the 
performance  of  various  DFT  functionals  was  studied  and 
compared  to  results  of  semiempirical  (INDO)  and 
Hartree  Fock  (HF)  methods.  (Isbom  et  al.,  2007)  In  that 
work,  molecular  structures  of  several  chromophores  were 
optimized  by  the  DFT  method  followed  by  a  calculation 
of  the  molecular  hyperpolarizability  using  semiempirical 
(INDO),  HF  and  DFT  methods,  the  latter  employing 
various  functionals  such  as  B3LYP  or  PBE.  The  authors 
conclude  that  the  relative  trends  of  hyperpolarizability 
from  various  computational  methods  are  similar,  and 
therefore  consistent  use  of  one  method  will  provide 
useful  guidance  for  an  experimental  program. 

The  visible  absorption  spectra  of  chromophores 
were  the  subjects  of  several  recent  studies  prompted  by 
the  importance  of  these  organic  molecules  in  the  dye 
industry.  Jacquemin  et  al.  (2008)  have  recently  reviewed 
the  performance  of  various  DFT  approaches  for 
predicting  A^x  in  the  visible  region  of  the 
electromagnetic  spectrum.  For  some  systems,  absolute 
errors  as  large  as  200  nm  were  reported  with  a  pure  DFT 
method  (PBE).  The  addition  of  25%  Hartree-Fock 
exchange  as  realized  in  the  PBE0-DFT  functional 
improved  results  considerably.  A  mean  absolute  error  of 
22  nm  and  a  maximum  deviation  of  less  than  100  nm 
were  reported  for  anthraquinone  and  azobenzene  dyes 
(Jacquemin  et  al.,  2008).  Dierksen  and  Grimme,  (2004) 
found  that  a  significantly  larger  amount  of  HF  exchange 
is  needed  to  improve  predictions  of  the  absorption 
spectra  for  systems  with  acetylenic  bonds. 

The  deficiencies  of  DFT  in  predicting  spectra  can  be 
associated  with  the  delocalization  error  of  approximate 
DFT  functionals  (Cohen  et  al.,  2008).  This  particularly 
affects  spectra  of  the  push-pull  molecules  investigated  in 
this  work.  This  error  can  be  partially  corrected  by 
including  an  HF  exchange  contribution,  either  constant 
as  stated  above,  or  asymptotically  approaching  full 
contribution  at  long  distances  between  electrons.  Long- 
range-corrected  DFT  (DFT-LC)  was  found  advantageous 
in  calculations  of  absorption  spectra  for  some  dye 
molecules.  (Jacquemin  et  al.,  2008) 


In  this  paper,  we  have  also  used  the  recently 
proposed  BNL  functional  (Livshits  and  Baer,  2007).  The 
BNL  functional  falls  into  this  DFT-LC  class,  and  it  is  of 
particular  interest  because  it  has  been  shown  to 
reproduce  exactly  the  charge-transfer  excitation  for 
model  systems.  Recently,  the  BNL  functional  was 
implemented  in  the  highly  parallel  NWChem  program 
(Andzelm  et  al.,  2008),  and  was  optimized  for  several 
types  of  chromophores. 

Finding  optimal  architectures  of  chromophores  that 
may  contain  numerous  substituents  in  various 
conformations  is  a  daunting  task.  The  number  of 
accessible  compounds  with  molecular  weight  below  850 
a.m.u.  (a  typical  size  of  EO  chromophores)  is  staggering 
at  about  1020°.  More  than  20  years  of  research  has 
established  many  rules  for  the  optimal  design  of 
chromophores,  effectively  reducing  our  design  space. 
Nevertheless,  surprising  new  architectures  of 
chromophores  are  appearing  at  increased  rates  with 
computational  chemistry  playing  an  important 
partnership  role  in  these  discoveries.  Efficient 
exploration  of  enormous  molecular  spaces  requires 
radically  new  optimization  schemes  that  search  directly 
for  the  property  optimum,  and  in  the  process  discover 
molecular  structure. 

The  typical  computational  search  algorithm  parallels 
the  experimental  paradigm,  by  first  generating  a  material 
structure  followed  by  a  calculation  of  its  properties.  With 
luck,  sometimes  the  results  are  satisfactory.  More  often 
than  not,  we  need  to  continue  repeating  the  search 
paradigm:  structure  —>  properties.  However,  if  this 
paradigm  could  be  inverted  into  properties  —> 
structure,  we  might  be  able  directly  to  discover  a  better 
material  with  the  required  properties.  In  this  work  we 
have  used  the  recently  proposed  inverse  molecular 
design  approach  (IMD).  (Keinan  et  al.,  2007)  The  IMD 
methods,  utilizing  the  concept  of  linear  combination  of 
atomic  potentials  (LCAP)  and  discrete  optimization 
approaches,  allows  for  the  highly-efficient  search  of 
optimal  properties  within  millions  of  possible 
chromophore  architectures. 

We  have  used  the  recently-developed  discrete 
optimization  branch-and-bound/tree-search  IMD 
algorithm  (Rinderspacher  et  al.,  2008)  to  search  for  the 
largest  hyperpolarizability  value  in  a  chemical  space  of 
chromophores.  In  this  paper,  we  have  extended  this 
algorithm  to  perform  a  constrained  IMD  search  (cIMD) 
that  allows  us  to  find  chromophores  with  large  values  of 
(3  and  simultaneously  visual  clarity.  The  constraint  was 
introduced  via  a  penalty  function  algorithm  with  a  soft 
penalty  weight. 

In  this  work,  we  analyze  the  properties  of  tolane 
molecules  and  develop  new  architectures  of  these  dyes 


that  enhance  their  P-value  while  retaining  visual 
transparency.  We  have  also  investigated  FTC-type 
chromophores  finding  substituents  that  enlarge  their 
visibility  window  without  affecting  significant  NLO 
properties  of  these  dyes. 


2.  METHOD 

The  chromophores  studied  in  this  paper  were 
calculated  using  either  DFT  or  the  semiempirical 
methods  AM1/ZINDO  with  Gaussian03,  NWChem  and 
CNDO  software.  The  structures  were  fully  optimized 
using  the  hybrid  DFT  functional  B3LYP  with  the  basis 
set  6-3 1G**  or  the  AMI  method.  The  DFT  UV-Vis 
spectra  were  predicted  from  singlet-singlet  excitations 
using  the  time-dependent  density  functional  theory 
(TDDFT)  method.  The  semiempirical  UV-Vis  spectra 
were  calculated  using  the  configuration  interaction 
method  with  single  excitations  (CIS)  within  the 
semiempirical  ZIND0/IND02  method.  The  solvation 
effects  were  considered  using  the  continuum  electrostatic 
models,  PCM  and  COSMO  as  implemented  in 
Gaussian03  and  NWChem,  respectively. 

The  EO  chromophores,  sought  in  this  work,  should 
possess  a  large  EO  coefficient,  which  is  a  function  of  the 
first  hyperpolarizability  ((3)  and  the  dipole  moment  (p). 
Therefore,  we  are  interested  in  maximizing  the  (3^  value, 
which  is  the  magnitude  of  the  irreducible  part  of  the  p 
tensor  projected  on  the  dipole  moment  of  the 
chromophore.  The  p  tensor  is  defined  as  a  nonlinear 
response  to  an  external  electric  field.  In  an  electric  field, 
F,  the  chromophore  energy  is 


E(F)  =  E(l-|iTF-l/2  aTF2-l/3!  (3TF3-... 


where  E0  is  the  energy  in  the  absence  of  the  electric 
field,  a  is  the  polarizability,  and  p  is  the  first 
hyperpolarizability.  The  hyperpolarizability  in  the 
direction  of  the  dipole  moment,  (3^  ,  is  then  calculated  as 

Pm-  =  (Zi  FiSj  Pijj)  /  M 

The  static  hyperpolarizability  p  was  calculated  using 
a  finite-field  method. 

2.1  LC-DFT  calculations 

The  BNL  functional  (Livshits  and  Baer,  2007)  based 
calculations  were  performed  using  the  NWChem 
computational  chemistry  package  developed  at  the 
Pacific  Northwest  National  Laboratory.  This  suite  of 
programs  has  been  designed  to  provide  the  capability  to 
run  large  scientific  molecular  simulations  on  massively 
parallel  and  scalable  computers.  In  the  BNL  approach, 
the  electron  repulsion  is  separated  into  long-  and  short- 
range  parts  and  given  as: 

—  =  —  erf  (yr)  +  —  erfc{yr) 
r 


where  y  is  an  attenuation  parameter.  The  short- 
range  part  is  treated  with  traditional  DFT  while  the  long- 
range  component  is  treated  with  exact  exchange.  It  can 
be  shown  that  for  large  y  the  theory  becomes  more  HF- 
like,  while  it  behaves  more  pure  DFT-like  for  small  y . 

We  have  recently  performed  an  optimization  of  this 
parameter  for  chromophores  containing  acceptors  with 
tricyano-motif  and  found  that  y  =  -0.35  leads  to  an 
accurate  prediction  of  absorption  spectra  with  a 
maximum  error  of  about  30  nm  for  the  most  intense 
absorption  peak,  Xmax.  The  BNL  implementation  in  the 
NWChem  program  results  in  a  highly  efficient  program 
to  calculate  absorption  UV-Vis  spectra  for  chromophores 
of  100  atoms  or  more  (Andzelm  et  al.,  2008) 

HF 
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Figure  2.  Fraction  of  HF  Exchange  vs.  inter-electron 
distance 

Figure  2  illustrates  the  role  of  exact-HF  exchange  in 
various  implementations  of  the  DFT  methodology.  The 
popular  density  or  gradient-density  dependent 
functionals,  LDA  or  PBE,  respectively,  do  not  include 
any  HF  exchange  and  therefore  exhibit  a  large 
delocalization  error.  The  constant  contribution  of  HF 
exchange  used  in  B3LYP  and  BHLYP  functionals  has  a 
localization  error  (Cohen  et  al.,  2008)  opposite  to  the 
delocalization  error  in  LDA.  Mixing  the  two  types  of 
functionals  causes  partial  cancellation  of  errors.  The 
DFT-LC  method  is  effective  for  this  error  cancellation  at 
large  electron-electron  distance  and  therefore  allowing 
for  accurate  calculations  of  absorption  spectra. 

Table  1  compares  the  performance  of  various  DFT 
methods  in  predicting  the  absorption  spectrum  of  dyes 
studied  in  this  paper.  For  each  dye,  two  main  absorption 
peaks  Vax  and  ^max-i  are  presented  together  with 
oscillator  strengths.  The  experimental  oscillator  strengths 
were  estimated  from  absorbance  data  and  normalized  to 
1.0  for  the  Xmax  peak.  The  6-3 1G*  basis  set  was  used  in 
DFT  and  HF  calculations  for  all  but  BHLYP,  where  the 
6-3 1+G*  basis  set  was  used. 

Adding  diffuse  functions  causes  a  red  shift  in  spectra 
by  up  to  20  nm.  To  compare  with  experimental  data,  the 
solvation  effects  of  chloroform  and  cyclohexanone  were 
included  for  the  tolane  and  FTC  dye,  respectively.  The 
solvent  causes  a  large  red  shift  of  the  >.max  peak  by  -90 
nm  for  the  FTC  dye  and  a  smaller  shift  of  -60  nm  for 
tolane  at  the  B3LYP/6-31G*  level  of  theory.  This  is  to 


r  r 


be  expected  because  the  dipole  moment  of  tolane  of  1 1.6 
Debye  is  much  smaller  than  that  of  FTC  at  22.7  Debye. 

Table  1.  UV-Vis  absorption  spectra  of  dyes  from  QM 
calculations.  Amax  and  'kmax-i  absorption  peaks  (in  nm)  and 
in  parentheses  oscillator  strengths  are  presented 


To \m  FTC 


Exptriritiil 

384(1.0) 

328(0.2) 

650(10) 

400(02) 

ZINDO 

407(1  j6) 

351(0.1) 

568(1.7) 

375(0.1) 

HF 

337(2.4) 

230(0.1) 

480(2.1) 

302(02) 

LDA 

965(0.4) 

520(0.9) 

811(13) 

553(0.7) 

PBE 

918(0.4) 

508(10) 

801(13) 

550(0.7) 

BHLYP 

409(20) 

312(06) 

573(20) 

368(0.1) 

83LYP 

586(0.7) 

404(1.4) 

690(1.7) 

464(05) 

M05 

520(00) 

38X14) 

657(1.7) 

440(0.4) 

BNL 

404(16) 

265(0.2) 

667(1.7) 

361(02) 

The  utility  of  the  computational  model  can  be 
judged  by  its  ability  to  reproduce  both  the  position  and 
strength  of  the  Amax  and  Amax_i  peaks  accurately,  thus 
providing  an  indication  of  the  width  and  absorption  level 
in  the  visible  part  of  the  spectrum.  Results  from  Table  1 
clearly  indicate  that  both  the  HF  and  the  DFT  methods, 
without  contributions  from  exact  exchange  (LDA,  PBE) 
are  useless  in  predicting  spectra  of  dyes.  The  DFT 
methods  with  some  HF-exchange  contributions 
(BHLYP,  B3LYP,  M05)  produce  erratic  results.  BHLYP 
is  apparently  useful  for  tolane  research,  while  M05 
predicts  well  the  Amax  peak  of  FTC.  The  DFT-LC 
methods,  such  as  BNL  implemented  by  us  (Andzelm  et 
al.,  2008)  seem  to  predict  consistently  the  main  spectral 
features  of  various  dyes. 


The  Amax  peak  represents  excitations  from  the  highest 
occupied  molecular  orbital  (HOMO)  to  the  lowest 
unoccupied  molecular  orbital  (LUMO).  This  is  a  charge- 
transfer  excitation  because  the  HOMO  is  localized  on  the 
donor  side  of  chromophores  while  the  LUMO  is  located 
on  the  acceptor  side.  The  significant  underestimate  of  the 
charge  transfer  HOMO-LUMO  gap  in  LDA  and  PBE, 


LUM0+1 


LUMO  ■ 


HOMO-1 
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Figure  3.  Valence  and  unoccupied  orbitals  of  tolane 


due  to  the  delocalization  error  (Cohen  et  al.,  2008)  is 
responsible  for  the  underestimate  of  the  excitation 
energy.  The  next  peak,  Amax-i,  is  due  to  n-  ti*  excitations, 
and  involves  orbitals  in  close  proximity  to  the  HOMO 
and  LUMO  levels.  Figure  3  shows  the  valence  and 
unoccupied  orbitals  of  the  tolane  molecule  that  dominate 
the  charge-transfer  and  7i-  7i*  excitations. 

2.2  cIMD 

The  inverse  molecular  design  method  optimizes 
molecular  property  by  minimization  of  its  gradient 
analogue.  In  our  case,  the  semiempirical  method  (AMI) 
is  used  to  obtain  geometries  and  then  by  changing 
substitutions  in  the  molecule,  the  approximate  gradient 
of  p  is  calculated  with  the  ZINOD/S  method.  Using  a 
discrete  analogue  of  the  steepest  descent  or  conjugated 
gradient  algorithm  one  can  find  the  optimal  structure  that 
maximizes  p.  The  key  idea  that  allows  for  “smoothness” 
of  the  property  surface  allowing  application  of 
continuous  optimization  techniques  is  the  definition  of 
the  electrostatic  potential  u(r)  of  a  molecule  as  a  linear 
combination  of  atomic  potentials  uR(r)  (LCAP): 

u(r)  =  Ik  0R  uV) 

The  0R  coefficients  are  continuous  variables  that 
define  the  admixture  of  an  atom  (or  group  of  atoms)  at 
position  R.  This  is  a  smooth  embedding  of  discrete 
Hamiltonians;  analogously,  the  properties  of  a  discrete 
set  of  molecules  can  be  embedded  in  a  smooth  space. 
Given  a  hierarchy  of  substitutions,  any  property  of 
interest  can  be  smoothly  interpolated.  Each  level  of  a 
hierarchy  of  substitutions  consists  of  a  list  of  substitution 
sites  and  a  set  of  subsequent  levels  for  each  site.  Then 
each  element  of  the  set  of  subsequent  levels  can  be 
identified  with  a  coefficient  from  0  to  1,  and  the  sum 
over  these  coefficients  for  each  set  must  equal  1 . 

In  the  following,  a  systematic  method  is  developed 
to  allow  the  application  of  continuous  optimization 
theory  to  discrete  problems.  This  is  achieved  by  an 
appropriate  embedding  of  the  discrete  set  of  property 
values  P  of  a  library  of  N  molecules.  Given  an 
enumeration  of  this  library,  where  each  molecule  is 
associated  with  a  number  s  and  its  property  value  Ps, 
log2N  bits  may  be  used  to  generate  a  continuous 
embedding.  For  example,  assuming  the  library  consists 
of  methane,  ethane,  propane,  and  butane,  in  just  that 
order,  then  it  is  possible  to  interpolate  between,  e.g.,  the 
ground  state  energies  (P=E)  of  these  4  molecules  using 
two  variables,  A0  and  A^  The  resultant  polynomial  is  just: 


E(Ao,  Ai)— E0(l-A0)(1  -Ai)h-E1A0(1  -Ai)+E2(1  -A0)Ai+  E3A0A1 

As  is  already  obvious  from  this  example,  the 
properties  of  the  polynomial  depend  on  the  ordering  of 
molecules.  Ordinary  interpolation  would  lead  to  a 


problematic  polynomial  of  degree  3.  Generalization  leads 
to  the  following  equations  including  derivatives,  where  s’ 
denotes  a  next  neighbor  of  s  at  bit  j : 

s=Z,2i-s(i),  s(i)£{0, 1 },  s'=s+(-l)s(i)-2s(i) 


performed  three  runs,  starting  from  three  different 
randomly  selected  structures.  All  the  runs  found  the 
expected  structure  as  shown  in  Fig.  4,  and  the  longest  run 
finishes  in  24  cycles,  well  below  the  theoretical 
maximum  of  200. 


P({Ai}i=1logN)  =  IS=0N  1  PsTTb=ilogN  ((1  -Ai)s(b|Ai1's,b)) 
A;=s(i),  dP/dAj=(-l)s<3)(Ps-Ps') 

This  highly  non-linear  polynomial  in  log2N  variables 
with  the  same  order  is  continuous  on  [0,  1]  log2N  and  it 
allows  the  adaptation  of  continuous  optimization 
methods  by  substituting  derivatives  via  simple 
differences.  Since  the  property  in  question  is  left  general, 
a  softly  weighted  penalty  function  was  used  to  enforce 
constraints.  The  following  algorithm  employed  in  this 
work  converges  due  to  the  local  convexity  of  P({Ai}: 

1.  Start  from  an  initial  molecule  and  mark  this 
molecule  as  current  best. 

2.  For  each  substitution  site  do  the  following: 

(a)  For  each  possible  substitution  at  the  current 
site,  attach  the  substitution  and  start  a  conformational 
search  and  compute  the  corresponding  p- values.  If  a 
constrained  optimization  is  performed  the  UV-Vis  is 
computed  and  a  penalty  is  applied  to  p.  The  result  of  the 
conformation  with  the  lowest  energy  is  recorded  for  this 
substitution. 

(b)  Compute  the  finite  difference  derivative  and 
choose  the  substitution  indicated  by  the  largest  change  in 
the  return  values  (p  for  unconstrained  optimizations)  for 
its  lowest  energy  conformation.  If  this  is  the  last 
substitution  site  proceed  with  3  else  continue  with  2  at 
the  next  site.  This  is  equivalent  to  choosing  the 
substitution  with  the  largest  p  for  its  lowest  energy 
conformation. 

3.  If  the  current  best  molecule  matches  the  last  best 
molecule  (structure,  energy,  property  value),  exit; 

otherwise  start  over  with  the  current  best  molecule  at 
step  2. 

Unlike  branch-and-bound  methods,  no  structures  are 
explicitly  excluded  from  the  search  space.  The  overall 
generality  of  the  algorithm  and  the  interpolation 
procedure  warrant  the  application  of  alternative 
optimization  methods  including  Monte  Carlo  or  genetic 
algorithms. 

This  method  was  tested  in  predicting  the  tolane 
structure  with  the  largest  p  value  as  a  function  of  four 
donor  (D)  substituents  and  two  acceptor  (A)  groups 
shown  in  Fig.  4.  The  best,  experimentally  known,  tolane 
architecture,  with  D=dimethyl-amine  and  A=nitro  group 
(Nguyen  et  al.,  1997),  is  also  displayed  in  Fig.  4.  In  this 
test  we  assume  that  the  donor  group  can  be  situated  in 
any  position  on  the  L  ring,  while  the  acceptor  group  can 
occupy  any  position  of  the  R  ring.  The  total  optimization 
space  is  then  5*4*5*2=200  molecules.  We  have 


HO  HS-  HJM-  (CH^N-  -NO,  -CM 

Figure  4.  The  optimal  structure  of  tolane  found  by  the 
IMD  method 


3.  RESULTS 


Conformations  of  chromophores  are  known 
significantly  to  affect  properties  such  as  the 
hyperpolarizability,  p.(Kinnibrugh  et  al.,  2006)  Much 
less  is  known  on  the  impact  of  conformations  on  the  UV- 
Vis  spectra  of  dyes.  One  can  expect  that  the  NLO 
properties  of  tolane  dyes  will  be  significantly  decreased 
due  to  almost  free  rotations  along  the  single  bonds 
connecting  phenyl  rings.  In  this  work,  we  will  analyze 
the  effect  of  tolane  bond  rotations,  and  propose  attaching 
substituents  in  the  meta-position  of  phenyl  rings  to 
increase  rotational  barrier.  We  will  initially  study  a  two- 
ring  tolane  (tolane-2),  p-amino-p’nitro  diphenylethyne. 


3.1  Properties  of  tolane-2 


The  structure  of  tolane-2,  rotational  barriers  and 
relative  changes  in  values  of  (3  and  as  a  function  of 
rotational  angle  along  the  single  bond  are  presented  in 


Figure  5.  Rotational  barrier  AE  (kcal/mol)  and  relative 
changes  in  Amax  and  (3  values  vs.  rotational  angle,  0 


Fig. 5.  The  maximum  barrier  at  90  deg  is  about  1.7  and 
1.5  kcal/mol  as  obtained  from  gas  and  solvent 
(chloroform)-phase  calculations  at  the  B3LYP/6-31G** 
level  of  the  DFT  theory.  Clearly,  any  rotamer  of  tolane-2 
can  easily  be  accessed  at  room  temperature.  The 
consequence  of  almost  free  rotation  is  a  decrease  in  the 
hyperpolarizability  (p)  value  by  a  factor  of  3.  The  effect 
of  free  rotation  on  the  UV-Vis  spectra  is  much  less 
pronounced.  The  ^max  peak  of  charge-transfer  character 
(384  nm)  dominates  the  spectrum  including  60  deg. 
rotamer.  Beginning  at  75  deg.  the  charge-transfer 
excitation  oscillator  strength  diminishes,  and  instead  the 
main  absorption  peak  has  a  n-n  character  at  a  much 
lower  absorption  of  300  nm.  Thus,  the  rotation  of  tolane 
dyes  contributes  to  the  significant  weakening  of  NLO 
properties,  and  also  results  in  a  blue-shift  of  the 
absorption  spectrum. 


Figure  6.  Dimer  of  tolane-2 

Tolane-2  forms  stable  dimer  structures  as  presented 
in  Figure  6.  The  B3LYP  level  of  theory,  so  successful  in 
predicting  structures  of  monomers,  fails  in  calculations 
for  dimmers,  because  of  the  lack  of  dispersion  forces  in 
the  DFT  functional.  The  MP2  method  yields  a  dimer 
structure  of  two  planar  tolane-2  monomers  at  a  distance 
of  about  3.68  A  measured  at  the  separation  between 
triple  bonds  on  the  monomers.  The  dimer  structure  is  due 
to  the  strong  electrostatic  interactions  of  two  opposite 
dipole  moments  of  11.3  Debye  each.  Due  to  opposite 
alignment  of  dipoles  the  total  moment  of  the  dimer 
cancels  out  and  the  hyperpolarizability  becomes 
negligible.  The  dimerization  of  tolane  molecules  has  no 
effect  on  the  UV-Vis  spectra. 

This  investigation  indicates  that  the  NLO  properties 
of  tolane-2  can  be  improved  if  the  free-rotation  of  phenyl 
groups  around  single  bonds  is  blocked  and  also  bulky 
substituents  are  attached  in  order  to  prevent  stacking  of 
tolane  monomers. 

3.2  Tolane-2  meta-substitutions 


Figure  7.  Hydrogen-bonding  interactions  between 
hydroxyl  and  aldehyde  groups  in  meta  positions  of 
tolane-2  for  rotamers  with  fixed  angle,  0  of  0,  30,  60  and 
90  deg;  AE  (in  kcal/mol),  ^max(in  nm)  and  P  (in  10"30  esu). 

Following  the  conclusions  of  the  preceding  section, 
we  have  attempted  to  constrain  free  rotations  around  the 
single  bonds  by  inserting  substituents  in  meta-positions 
of  phenyl  rings  that  strongly  interact  and  raise  the  barrier 
for  rotation.  Figure  7  shows  an  example  of  a  modified 
tolane-2  structure  with  high  rotational  barrier  due  to 
strong  hydrogen  bonding  interaction  between  the 
hydroxyl  and  aldehyde  groups.  The  barrier  for  rotation 
(AE)  is  18.7  kcal/mol.  For  this  rotamer,  the  main  Lmax 
peak  of  301  nm  has  n-n  character,  and  is  significantly 
blue  shifted  in  comparison  to  the  lowest  energy  rotamer. 
As  expected,  the  hyperpolarizability  is  significantly 
decreased  due  to  the  rotation.  The  pair  of  substituents  - 
hydroxyl-aldehyde  leads  to  an  excellent  rotational 
barrier;  however,  the  main  absorption  peak  of  about  450 
nm  makes  this  structure  too  red  for  our  application. 


We  have  examined  numerous  pairs  of  substituents 
searching  for  structures  with  low  absorption  not 
exceeding  430  nm  and  with  strong  hydrogen  bonding 


L  R 


AE 


A 


max 


-ch2nh2 

-ch2oh 

-cncoch3 

-COOCH3 

-nh2 

-CHOCH3 

-nh2 

-CH50H 

-nh2 

-CN 

-nh2 

-cf3 

“OH 

-COH 

-NH2 

-COOH 

-nh2 

-COCH3 

0  5  10  15  20  400  430  400 


Figure  8.  Rotational  barriers  (AE)  and  for  tolane-2 
with  substituents  in  meta-positions 


stabilizing  internal  rotations.  The  results  are  collected  in 
Figure  8.  The  amino  groups  (-NH2)  are  strongly  bonded 
to  ketone  (-COCH3),  acid  (-COOH)  and  alcohol  (- 
CH2OH)  substituent  groups.  The  hydroxyl  group  (-OH) 
interaction  with  aldehyde  (-COH)  and  also  amine 
interaction  with  cyano  groups  (-CN)  also  yield 
significant  rotational  barriers.  Unfortunately,  for  all  these 
cases  the  main  absorption  peak  is  well  above  430  nm 
indicating  an  opaque  chromophore  that  is  unacceptable 
for  our  application.  The  most  transparent  chromophores 
with  methylamine-,  methoxy-,  or  amino-ether 
substituents  have  rather  low  rotational  barrier.  An  amido 
(-NHCOCH3)  -ester  (-COOCH3)  pair  was  found  to 
produce  a  rather  large  barrier  to  rotation  at  7.0  kcal/mol 
and  an  acceptable  position  of  the  ^max  peak  at  424  nm. 

Figure  8  displays  the  lowest  energy  conformers  of 
amine-ketone,  and  also  the  amido-ester  structures.  The 
rotational  barrier  can  be  improved  by  attaching 
substituents  that  form  multiple  hydrogen  bonds  resulting 
in  structures  with  excellent  visual  properties  (^max=  -400 
nm)  and  significant  barriers  to  rotation  of -17  kcal/mol. 

3.3  Tolane-3  meta-substitutions 


Figure  9.  Tolane-3  with  amido-ester  motif  forming 
multiple  hydrogen  bonds 

(3  =  225  10"30  esu.  The  double  hydrogen  bonds  between 
side  chains  of  adjacent  phenylene  units,  stabilize  the 
internal  rotation  yielding  a  barrier  of  10  kcal/mol.  We  are 
currently  working  on  a  synthesis  of  amido-ester 
functionalized  tolane  dyes  with  multiple  internal 
hydrogen  bonds.  The  synthesis  of  p- 
phenyleneethynylenes  with  internal  hydrogen  bond  has 
been  reported.  (Hu  et  al.,  2008)  However,  these  are  not 
push-pull  type  EO  molecules  of  interest  in  this  work 
(Cross  and  Davis,  2008),  and  the  synthesis  of  dyes  with 
the  amido-ester  motif  as  shown  in  Fig.  9  is  considerably 
more  challenging. 


We  have  used  cIMD  method  to  explore  the  best 
substitution  sites  on  tolane-3  that  maximize  p  values  and 
maintain  transparency.  The  optimization  space  consisted 
of  several  substitution  groups  located  in  A-D’ sites  shown 
in  Fig.  4.  Our  first  run  (I)  included  5  substitutions  (- 
CH2NH2,  -CH2OCH3,  -COOCH3,  -NHCOCH3,  -H) 
visiting  all  sites  A-D  for  a  total  of  58  =  390,625 
configurations.  The  second  run  (II)  used  a  limited 
number  of  substitution  sites  (A,  B,  C,  D)  for  a  total  of 
54=625  configurations.  The  third  run  (III)  used  9 
additional  substituents  (-NH2,  -CH3,  -OH,  -COCH3,  - 
OCH3,  -CF3,  -COH,  COOH,  -CH2OH)  located  at  the  A, 
B,  C,  D  sites  creating  an  optimization  space  of  144  = 
38,416.  The  analysis  of  results  was  limited  to  the  20 
conformations  with  largest  value  of  P  that  also  form  at 
least  two  hydrogen  bonds  stabilizing  rotations  around  the 
acetylene-aryl  single  bonds  and  have  ^ax  <  420  nm. 

We  find  that  amido-ester  pairs  occur  more  often  in 
this  set  of  preferable  configurations.  A  small  populations 
of  methyloamino-ester  and  amido-ether  pairs  was  also 
found.  Amide  groups  were  mainly  located  in  meta 
positions  A,  A’,  C  and  C\  Test  number  II  confirmed  that 
amide-ester  pairs  at  A-B  and  C-D  provide  the  best  p 
value.  Optimization  III  led  to  various  amine  complexes 
with  CF3  groups  that  could  not  be  accepted  because  of 
the  red-shifted  spectra. 

Inspired  by  the  cIMD  simulations,  we  have  built  a 
system  based  on  an  amido-ester  motif  (see  Fig  9)  for 
tolane-3  in  A-B  and  C-D  positions  (Fig.  4)  that  has 
excellent  transparency  ^max  =414  nm  and  a  high  value  of 


3.4  FTC  substitutions 

The  FTC  chromophore  (Fig  1)  has  several 
conformers  that  are  close  in  energy;  however,  their 
hyperpolarizabilities  vary  by  as  much  as  40%.  FTC  could 
still  be  an  attractive  chromophore  for  Army  applications 
if  not  for  its  poor  transparency.  The  ^max_value  is  around 
650  nm  blocking  effectively  the  yellow  visibility  region 
of  the  spectra.  In  this  work  we  are  concerned  with 
designing  an  FTC  architecture  that  red-shifts  the  7^^  - 
peak. 

The  energies  of  optimized  structures,  for  the  lowest 
six  FTC  conformers  are  close,  within  2  kcal/mol.  We 
also  confirm  that  the  P  values  may  decrease  significantly 
(from  305  to  193  xlO"30esu),  (Kinnibrugh  et  al.,  2007); 
however  we  find  that  this  is  not  the  case  for  the  Xmax  - 
peak.  The  variations  in  the  value  of  the  main  absorption 
peaks  are  within  1 5  nm  of  the  7^^  for  the  most  stable 
conformer  We  conclude  that  conformation  changes  of 
FTC  only  slightly  affect  the  UV-Vis  absorption  spectra. 


Figure  10.  The  ethoxy- substituted  FTC  dye  with 
improved  ^max 


We  performed  a  cIMD  search  of  the  best 
substitutions  in  the  lowest  energy  FTC  conformer  that 
maximizes  the  (3-value  and  constrains  ^ax  to  be  higher 
than  the  original  value.  The  search  space  includes  four 
substitution  sites  (marked  with  X=A,  B,  C,  D  in  Fig  10), 
choice  of  donor,  and  exchanging  one  of  the  methyl 
groups  of  an  acceptor  with  CF3.  The  following  groups 
were  considered  for  the  X-sites:  hydroxyl,  methyl, 
methoxy,  methylthioether,  amide,  ester  and  hydrogen. 
We  have  also  considered  amine,  methyl-amine  and  ethyl- 
amine  as  donors.  The  total  optimization  space  contained 
14,406  molecules  and  our  search  completed  in  62  cycles. 
Analysis  of  the  top  12  structures  revealed  that  the 
preferable  architecture  has  a  CF3  group  at  the  acceptor 
site  and  ethyl-amine  as  a  donor.  The  best  structure  shows 
methoxy  substitutions  at  the  A,  C  and  D  sites.  The  A  site 
is  predominantly  populated  with  a  methoxy  group  (one 
occurrence  of  hydroxyl  and  amide  groups  was  found). 
The  geometries  of  these  structures  were  then  refined  with 
DFT  and  calculations  of  UV-Vis  and  hyperpolarizability 
confirmed  the  findings  of  the  IMD  method.  The  best 
structure  is  presented  in  Figure  10.  A  significant  red-shift 
of  the  ^max-peak  of  about  15%  (to  740  nm),  and  even  a 
slight  increase  in  (3-value  (327x1 0"30  esu)  was  found. 

CONCLUSIONS 

In  this  work,  we  have  presented  a  new  constrained 
inverse  molecular  design  algorithm  (cIMD)  and  the 
long-range  corrected  density  functional  method  (LC- 
DFT).  We  have  shown  that  these  methods  facilitate  the 
discovery  of  new  chromophores  for  Army  applications. 
Specifically,  we  have  found  that  the  self-assembly  of 
hydrogen-bonded  substituents  in  the  meta  position  of 
tolanes  prohibits  internal  rotations,  thus  improving 
significantly  the  (3-value.  The  amido-ester  motif 
involving  multiple  intramolecular  hydrogen-bonding  was 
found  as  the  best  compromise  between  requirement  of 
achieving  high  rotational  barriers  and  low  UV-Vis 
absorption.  The  FTC  chromophores  were  also  modified 
by  including  CF3  and  methoxy  substituents  that  enlarge 
the  visibility  window  without  affecting  the  high  (3- values. 
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